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MAGNETIZED ASTROPHYSICAL OUTFLOWS: 
CRADLE TO GRAVE, SOURCE TO EFFECT 

A. Frank, T. Gardiner 1 and T. Lery 2 
RESUMEN 

El resumen sera traducido al espanol por los editores. 

We discuss the propagation of radiative MHD jets and outflows focusing on outflows driven by magneto- 
centrifugal rotators. Our goal is to link the properties of the jets with the physics of the sources which produce 
them. We find that density and magnetic field stratification (with radius) in jets from magnetized rotators 
leads to new behavior including the development of a dense inner jet core and a low density collar. We also 
report on more general studies of ambipolar diffusion and field geometry in pulsed jets. Finally we describe 
a new work designed to study the effects of magnetized winds on circumstellar environments appropriate to 
YSOs and PNe. 
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We discuss the propagation of radiative MHD jets and outflows focusing on outflows driven by magneto- 
centrifugal rotators. Our goal is to link the properties of the jets with the physics of the sources which produce 
them. We find that density and magnetic field stratification (with radius) in jets from magnetized rotators 
leads to new behavior including the development of a dense inner jet core and a low density collar. We also 
report on more general studies of ambipolar diffusion and field geometry in pulsed jets. Finally we describe 
a new work designed to study the effects of magnetized winds on circumstellar environments appropriate to 
YSOs and PNe. 
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1. INTRODUCTION 

Hypersonic collimated flows (jets) are a ubiq- 
uitous phenomena in astrophysics occurring in 
wide variety of environments including young stars 
(YSOs) (Reipurth 1997) and evolved stars (PNe) 
(Soker & Livio 1994). In almost all cases magnetic 
fields are expected to play a key role in launching and 
collimating these outflows. In particular, magneto- 
centrifugal processes associated with magnetized, ro- 
tating stars (Bogovalov & Tsinganos 1999) and/or 
accretion disks (Pudritz 1991) are believed to lift 
material out of the gravitational well and provide 
at least some confinement, shaping the wind into a 
jet. There has been considerable effort in the study 
of collimated outflows over the last two decades both 
in terms of jet simulations and the study of magneto- 
centrifugal launching. Unfortunately, direct links 
between the observable jets (scales of > 10 15 cm) 
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and the collimation processes (scales of < 10 13 cm) 
have yet to be established. Indeed, until recently the 
majority of jet propagation simulations have been 
purely hydrodynamic. Thus the role magnetic fields 
play in establishing the properties and behavior of 
radiative jets (appropriate to YSOs and PNe) re- 
mains an open issue. In addition, establishing links 
between the near-field (close to the star) magneto- 
centrifugal processes and the far-field jet behavior 
which can be readily observed must also be estab- 
lished. 

In this contribution we describe the MHD be- 
havior in radiative jets and outflows. The goal of 
our studies is to articulate how magnetic fields effect 
radiative jet/outflow behavior and to provide links 
between magneto-centrifugal processes occurring at 
the source with observable properties of the jets and 
outflows. 

2. JET STRUCTURE AND THE PROBLEM OF 
INITIAL CONDITIONS 

Magnetic fields imposed on the flow during the 
launching process may bleed out the beam by am- 
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bipolar diffusion (where ions and the field slip past 
neutrals. Frank et al (1999) have shown that the 
ambipolar diffusion timescalc in YSO jets can be 
written as 

where (3 is the ratio of gas to magnetic pressure, 
Q{fi) — (3/(/3 + 1), and rij refers to the neutral frac- 
tion in the jet beam. Since the largest scale YSO 
jets have dynamical times of order a few 10 4 to 10 5 y 
(Reipurth 1997) our result shows that fields will re- 
main in the beam for much of the jet's lifetime. The 
characteristic ambipolar diffusion timescale of 10 4 y 
is also suggestive. It's approximate equivalence with 
the age of parsec-scale jets may indicate that these 
structures remain intact as long as the collimating 
fields remain in the beam. This issue needs further 
study. 

Establishing initial equilibria for MHD jet simu- 
lations is non-trivial. This is not the case for hy- 
drodynamic jets where the required force balance 
across the jet and ambient medium interface allows 
for the use of so-called top-hat profiles (i.e. the hydro 
variables are constant across the jet cross section). 
Such distributions may not be tenable in MHD jet 
studies. The difficulty can be seen by decomposing 
the Lorentz force into a tension term and a pressure 
term. Fi oc -VB 2 + 2(B • V)B. In a steady, cylin- 
drically collimated jet only and B z components 
of the field are possible. Jets with purely longitu- 
dinal fields, B = B z , can be easily set in pressure 
balance with the environment and top-hat profiles 
may be used. Toroidal or helical field geometries re- 
quire more complicated initial conditions unless the 
field is assumed to take on a force free configuration. 
If the field is not force free, MHD jets must have 
variable distributions of gas pressure and, perhaps, 
other variables in order to balance the hoop stresses 
associated with the tension force. 

Faced with the problem of initial conditions re- 
searchers studying radiative MHD jets have, in gen- 
eral adopted one of two strategies: (i) use force-free 
fields (Cerqueira et al 1998; Cerqueira et al 1999); 
(ii) use ad-hoc gas pressure and magnetic field dis- 
tributions configured to be in initial force balance 
(Frank et al 1999; Stone & Hardee 2000; O'Sullivan 
& Ray, 2000). The results of these studies for both 
steady and time- variable (pulsing) jets reveal a num- 
ber generic features. 

• Jets with purely longitudinal geometries do 
not show propagation characteristics which dif- 
fer significantly from the hydrodynamic case 



(Cerqueira et al 1998; Gardiner et al 2000). 
The presence of poloidal fields does however al- 
low for the possibility of field reversals at the 
head of the jet which are likely to be unstable 
to reconnection (Gardiner et al 2000). The 
transfer of magnetic energy into thermal energy 
at reconnection sites may alter the emitted spec- 
trum away from that associated with shocks. 
Given the importance of shock diagnostics for 
interpretation of spectra the additional source 
of excitation provided by reconnection requires 
further study. 

• Jets with a toroidal field component will be sub- 
ject to strong hoop stresses especially in the re- 
gion between the jet and bow shocks (Clarke 
et al 1986; Frank et al 1998). In 2.5-D ax- 
isymmctric calculations pinch forces associated 
with the toroidal fields lead to the development 
of streamlined, high field nose-cones. Recent 
simulations have shown such nose-cone may be 
unstable in 3-D (Dal Pino, these proceedings). 
While this result still needs further study, it may 
be correct for jets with top-hat density profiles. 
As we will show in the next section, more re- 
alistic initial conditions for jets from magneto- 
centrifugal rotators may not have unstable nose- 
cones. 

While the studies discussed above have yielded 
significant progress in our understanding of heavy, 
radiative MHD jet behavior, the initial conditions 
arc still unconnected to the processes believed to 
create the jets. We turn to this issue in the next 
section. 

3. JETS FROM KEPLERIAN ROTATORS 
Recently we have simulated jets whose initial 
conditions jets are taken directly from a (sim- 
plified) model of the magneto-centrifugal launch- 
ing/collimation process. The model, known as the 
Given Geometry Method (GGM) allows asymptotic 
MHD jet equilibria to be linked to the properties of 
a rotating source (Lery et al 1999; Lery & Frank 
2000). The GGM assumes a time- independent, ax- 
isymmetric flow. It further simplifies the problem 
of magneto-centrifugal launching/collimation by as- 
suming that the nested magnetic flux surfaces defin- 
ing the flow possess a shape which is known a priori 
inside the fast critical surface (the locus of points 
beyond which the flow is kinetic energy dominated). 
The flux surfaces are assumed to be conical and, as 
an additional simplification, an equilibrium across 
the surfaces is assumed at the Alfven point. The 
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final asymptotic solution for the collimated jet is 
solved by assuming pressure balance with the am- 
bient medium. 

In terms of dynamics, The most important prop- 
erties of the jets derived via the GGM are the radial 
variations in the density p(r) and the toroidal com- 
ponent of the magnetic field B^r). In particular 
when the source is composed of a rapidly rotating 
disk truncated some distance from a rigidly rotating 
star, the emitted jets can have strong density strat- 
ifications, ie. a high density axial core surrounded 
by a lower density collar. The strongest toroidal 
field is located at the boundary of the core and col- 
lar creating a magnetically confined jet-within-a-jet 
structure. Note that the bulk of the jet's momentum 
resides in the core. Hence we expect this portion of 
the beam to penetrate more easily into the ambient 
medium during the jet's propagation while the collar 
will be more strongly decelerated. 

Figure 1 (top) shows a frame from a jet simula- 
tion with a high core to collar density ratio (~ 100, 
see Frank et al (2000)). As expected the jet core 
propagates faster through the grid then the outer 
collar. In these simulations, (which do not include 
radiative losses), the collar is diverted by shocks at 
the initial head of the jet. The collar peels away as 
the jet propagates down the grid. A distinct nose- 
cone forms at the head of the jet core due to its 
strong confining toroidal magnetic field. Nose-cones 
associated with this high density core may not be un- 
stable. Also of particular interest are the apparent 
pinch mode (m = 0) instabilities seen in the core. A 
stability analysis of the initial conditions effectively 
predicts the wavelength of these modes (A = .3Rj) 
and suggests these they are current driven rather hy- 
drodynamic Kelvin-Hclmholtz instabilities. 

Since many jets are seen to be episodic we have 
can also use the GGM method to derived a sequence 
of jet solutions appropriate to a time variable flow. 
Figure 1 (bottom) shows the density in a radiative 
pulsed jet simulation (Gardiner, Frank & Lery 2001) 
from a Keplerian rotator. Here again we see the 
core/collar structure which now persists in the inter- 
nal working surfaces. The collar is not seen to peel 
away from the core in these models, a result which 
may be due to the pulsing or to the lower core to 
collar density ratio obtained from the GGM model 
in this case. The eventual goal of these studies will 
be to model jets emitted by sources undergoing some 
form of FU Ori type outburst and hence connect the 
properties of variable jets to theoretical models of 
variations at the source. 
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Fig. 1. Top: Grey-scale map of the density in steady 
jet calculation from a disk+star rotator (see Frank et 
al 2000 for details). Bottom: Pulsed jet simulation. 
Initial conditions derived from a sequence of GGM Kep- 
lerian rotator models (Gardiner, Frank & Lery 2001). 

4. THE MAGNETIC GEOMETRY OF PULSED 
JETS 

The most general magnetic geometry for cylindri- 
cally collimated jets is a helical field. Indeed, this is 
the geometry which is to be expected from magneto- 
centrifugal launching. Since many jets are observed 
to episodic, it is of interest to understand how vari- 
ations in the flow at the source will effect the sub- 
sequent downstream field geometry of the jets. In 
Gardiner & Frank (2001a) a simple, universal re- 
lation was derived for the field configuration which 
would evolve in a helically magnetized, pulsed jet. 

Consideration of the governing flow equations re- 
veals that the toroidal field component will evolve 
in manner similar to the gas density. Thus, down- 
stream compressions produced by internal working 
surfaces will also lead to compression in the toroidal 
field. Likewise, regions of gas rarefactions will also 
become regions of weak toroidal field. Gardiner & 
Frank (2001a) derived an expression for the ratio of 
the toroidal to poloidal field components at a time t 
based on a Burger's equation analysis, 

Bff, _ B^fl / 1 \ 

s7 ~ "bT \i-k(t-t )) 

where t a is the time when the gas parcel was ejected 
and k is the ratio of the derivative of the jet ejection 
velocity at t Q to the value of the velocity at that time 

* = VjltoY/Vjito). 

The relation above shows that pulsed helical 
jets will develop alternating regions of toroidal and 
poloidal domination in the beam. High density knots 
associated with the internal working surfaces will be 
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Fig. 2. Logio(B<p/B<f,o) at two different times in a pulsed 
radiative MHD jet simulation with an initial helical field. 
For details see Gardiner & Frank (2001a). 

toroidally dominated and the inter-knot rarefied re- 
gions will be poloidally dominated. Figure 2 shows 
a plot of B$ at two different times from a pulsed 
jet simulation. The development of strong field re- 
gions is clearly seen as the pulses steepen into inter- 
nal shocks. 

5. MAGNETIZED WIND BLOWN 
BUBBLES:YSO AND PNE 

In previous sections we discussed the behavior of 
fully collimated jets. In both YSOs and PNe, how- 
ever, the outflows may include a significant wide- 
angle magnetized wind component. Numerous stud- 
ies have shown that the pure hydrodynamical inter- 
action of a wide-angle wind with an asphcrical ambi- 
ent density distribution can lead to highly collimated 
bipolar flows (Icke et al 1992; Dclamarter, Frank & 
Hartmann 2000). Since the environment surround- 
ing both YSOs and PNe is expected to include an 
aspherically shaped density distribution, the addi- 
tional effect of magnetic fields in shaping the out- 
flows is of great interest. 

A number of authors have explored a model 
whereby a weak toroidally magnetized wind is col- 
limated by hoop stresses after passage through a 
wind shock (Chevalier & Luo 1994; Garcia-Segura 
et al 1999). This so-called magnetized wind 
bubble (MWB) model shows some promise in that 
highly collimated outflows can be obtained even 
when purely hydrodynamic flows would only yield 
rather wide bipolar lobes. Precession of the sources' 
magnetic axis also allows for point-symmetry, (ob- 
served in many PNe), to be imposed on the flows 
(Sahai 2000). 

Recently however (Gardiner & Frank 2001b) 



Fig. 3. Simulation of MWB model in a collapsing en- 
vironment appropriate to YSOs. The simulation be- 
gins with a wind with properties similar to those from 
magneto-centrifugal flows. 

have shown that the MWB fails to account for colli- 
mation which would occur before the wind is shocked. 
The toroidal field in the wind produces an unbal- 
anced collimating force which can redirect the wind 
as it freely expands. Thus a proper treatment of this 
model must include the full history of the wind start- 
ing close to the source. It is also noteworthy that the 
MWB model is consistent with a wind accelerated by 
some means other than magneto-centrifugal launch- 
ing. 

Figure 3 shows a simulation of a magnetized wind 
blown bubble appropriate to YSOs (Gardiner, Frank 
& Hartmann 2001). In this model a pre- collimated 
wind with a pole to equator density ratio of 10 is 
driven into a collapsing sheet-like environment. The 
initial conditions capture many of the salient proper- 
ties of magneto-centrifugal flows such as the X-wind 
(Shu et al 1994; Matzner & McKee 1999). The fig- 
ure shows a strongly collimated flow, including a jet, 
forming along the axis. Note that though the flow 
begins with a strong axial condensation it is not visi- 
ble until after the wind passes through the shock. As 
with previous studies of the MWB model this occurs 
due to strengthening of hoop stresses associated with 
the toroidal field after the wind has been shocked. 

6. CONCLUSIONS 

Magnetic fields offer a theoretically attractive 
means for producing and shaping much of the out- 
flow behavior observed in both young and evolved 
stars. The obvious failing of magnetic models rests in 
the unfortunate reality that the fields themselves are, 
in most cases, unobservable with current techniques. 
Thus the presence and effect of the these important 
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players on the stage of stellar evolution must be in- 
ferred by other means. Currently it is not clear if 
YSO jets, the structures most likely associated with 
stellar or disk magnetic fields, show any behaviors, 
such as nose-cones, which can be clearly linked to 
the presence of magnetic stresses. It would be an 
ironic situation indeed if the most favored model of 
jet launching and collimation leaves us with no sig- 
nature with which to confirm its veracity. 

It is still too early to know if simulation stud- 
ies can provide proxy links between outflows and 
magnetic fields. We believe, however, that the use 
of initial conditions which come directly from MHD 
launching and collimation models should provide the 
quickest route to a final answer. 
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